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Abstract. We present a catalogue of astrophysical data of 520 Galactic open clusters. These are the clusters, for which at least 
three most probable members (18 on average) could be identified in the ASCC-2.5, a catalogue of stars based on the Tycho-2 
observations from the Hipparcos mission. We applied homogeneous methods and algorithms to determine angular sizes of 
cluster cores and coronae, heliocentric distances, mean proper motions, mean radial velocities and ages. For the first time 
we derive distances for 200 clusters, radial velocities for 94 clusters and determine ages of 196 clusters. This homogeneous 
new parameter set is compared with earlier determinations. In particular, we find that the angular sizes were systematically 
underestimated in the literature. 
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1. Introduction 

Although open clusters are typical representatives of the 
Galactic disk population, a systematic investigation of their na- 
ture, size, number of members and age is hampered by the in- 
homogeneity of the data. On one hand, general bibliographic 
catalogues derived from reviews of the literature are available. 
Examples of this category are the catalogues by Lynga ( 1987 1, 
Ruprecht et al.( U981> . Dias et al. (20041, and the data base 
WEBDA 1 of Mermilliod. In these collections general parame- 
ters of the clusters are taken from the literature as published. 
These collections are indispensable sources for work in this 
field. However, their parameters cannot be used for system- 
atic studies or intercomparisons between clusters. On the other 
hand, there are uniform lists of a few hundreds of clusters, 
the parameters of which were derived either from coherent 
photometric studies (e.g. Becker and Fenkart fl971l Janes and 
Adrer [T9"g2l Loktin et al. EOTFTl Tadross UMH . or from their 
kinematics (e.g. Baumgardt et al. l2000l Dias et al. l200!l 120021 . 

A few years ago we started the project to determine mem- 
bership and principal parameters of open clusters by use of ob- 



Send offprint requests to: R.-D. Scholz 
1 http://obswww.unige.ch/webda/ 



jective methods and algorithms. In order to get such an un- 
biased view onto the system of open clusters of our Galaxy, 
one has to start with an homogeneous sky survey obtained with 
no particular emphasis on open clusters. One result of ESA's 
Hipparcos mission, the Tycho-2 catalogue (H0g et al. 2000 1, 
fulfills this necessary criterion, but its contents alone were not 
sufficient for our purpose. 

So, as a first step, a catalogue of 2.5 million stars with 
proper motions in the Hipparcos system and B, V magnitudes 
in the Johnson photometric system was compiled and supple- 
mented with spectral types and radial velocities if available. 
The resulting catalogue, the All-Sky Compiled Catalogue of 
2.5 Million Stars (ASCC-2.5, Kharchenko I200H can be re- 
trieved from the CDS 2 , a detailed description of the catalogue 
can be found in that paper or in the corresponding ReadMe file 
at the CDS. 

In a second step, we used the ASCC-2.5 to identify 
known open clusters and compact associations listed in the 
Lund Catalogue of Open Clusters (Lynga 1987), the Optically 
Visible Open Clusters Catalog (Dias et al. 2004 referred 
hereafter as DLAM), and the Catalogue of Star Clusters and 
Associations (Ruprecht et al. U981> . Applying an iterative pro- 
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cedure or cluster membership determination based on proper 
motion, photometric and spatial criteria to the ASCC-2.5 data, 
we could identify 520 of about 1700 known clusters. After re- 
estimating the positions of cluster centres and cluster sizes, 
about 150000 ASCC-2.5 stars were selected in these 520 clus- 
ter areas, and membership probabilities were computed for 
each star in this list which we called the Catalogue of Stars in 
Open Cluster Areas (CSOCA). The procedure of cluster mem- 
bership determination is briefly described below, a detailed in- 
formation is given in Kharchenko et al. ( 2004b I, referred here- 
after as Paper I. CSOCA is a catalogue of stars which now gives 
us the basis for deriving uniform structural (location, size), 
kinematic (proper motions and spatial velocities) and evolu- 
tionary (age) parameters for open clusters in the wider neigh- 
bourhood of the Sun. The results are published in the Catalogue 
of Open Cluster Data (COCD) and the Open Cluster Diagrams 
Atlas (OCDA). 

The present paper is the third step in our long-term project, 
and it is structured as follows: in Sec.[2]we briefly describe the 
basic data and the procedure of member selection we devel- 
oped. In Sec. 13 the parameters describing the spatial-structure 
of the clusters are discussed. Sec.0]is devoted to the kinemat- 
ics of the clusters. In Sec.|5]we describe the method applied for 
deriving ages of the clusters. Concluding remarks are given in 
Sec.|6] In the appendices we describe the format of COCD and 
OCDA. 

2. The cluster list, and cluster membership 

The input data for this study come from the ASCC-2.5 and 
cover 513 open clusters and 7 compact associations (i.e. the 
associations with a sufficiently high projected density of stars, 
so that the developed general procedure of member selection 
works properly). Due to a relatively bright limiting magnitude 
(V ~ 14 mag) of the ASCC-2.5, our sample does not include 
faint (and generally remote or highly obscured) open clusters. 
Also, the two nearest extended clusters, the Hyades and Cr 285 
(the UMa cluster), are missing in our list since they require a 
special and more sophisticated technique of membership deter- 
mination. However, from comparison with DLAM, our sample 
is sufficiently complete for clusters up to 1 kpc and can be used 
for the study of cluster parameters and properties of the local 
population of open clusters. 

For each cluster, the membership determination is based on 
a comprehensive common analysis of several diagrams derived 
with ASCC-2.5 data: a sky chart referring to the cluster, the 
radial distribution F{r) of the projected stellar density, a vec- 
tor point diagram (VPD) of the proper motions, the magnitude 
dependence of the proper motion components, and a colour- 
magnitude diagram (CMD). 

As a first approximation of the member selection algorithm 
applied, we consider all ASCC-2.5 stars within an area of lxl 
sq. deg around the cluster centre with coordinates taken from 
DLAM. Using proper motion and photometric criteria, we then 
separate field stars and cluster members and adopt the point 
of maximum surface density of cluster members as the new 
approximation of the location of the cluster centre. From the 
analysis of the surface cluster density F(r), a projected radius 



Fig. 1. Illustration of the member-selection algorithm for the 
open cluster Stock 2. Arrows show the sequence of the basic se- 
lection phases presented by rectangles. Panel (a) is a sky map 
of the cluster neighbourhood. Panel (b) shows radial profiles 
of the projected density, panel (c) is the colour-magnitude di- 
agram, panel (d) is vector point diagram, and panels (e) and 
(f) are "magnitude equations" (ji x>y - V relations; they are use- 
ful as a proper motion check since proper motions of cluster 
members should not depend on their magnitudes). Additional 
explanations for each panel are given in AppendixlBl 

of the cluster is then derived. The resulting spatial parameters 
are used for the next steps in the iteration. The iterations are 
stopped when the cluster member list does no longer change. 
As a rule, two iterations are sufficient to establish cluster mem- 
bership and determine the structural and kinematic parameters. 
However, some clusters with sparse structure or initially er- 
roneous central coordinates require three and more iterations. 
The member selection pipeline is illustrated in Fig^fd" the 
case of the open cluster Stock 2 (see Paper I for more details). 

For stars, the probability of belonging to a cluster is calcu- 
lated as a measure of a deviation either from the cluster mean 
proper motion (kinematical probability), or from the Main 
Sequence (MS) edges (photometric probability). Stars deviat- 
ing from the reference values by less than one cr rms are classi- 
fied as most probable cluster members (lcr-members, i.e., with 
a membership probability P > 61%). Those falling in semi in- 
tervals [l<x,2cr) or [2<x,3cr) are considered as possible members 
(P = 14 - 61%) or possible field stars (P = 1 - 14%), re- 
spectively. Stars with deviations larger than 3<x are regarded as 
definite field stars (P < 1%). As a rule, all cluster parameters 
were determined by use of data on the most probable cluster 
members. 

3. Spatial structure parameters of the clusters 

3.1. Distances and extinction 

Although about 9 700 of the CSOCA stars with cluster 
membership probabilities higher than 14% (i.e. lcr- or 
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Fig. 2. Comparison of distances d and colour-excesses 
E(B - V) derived by Loktin et al. (EtUll from 3-colour 

photometry with values determined on the basis of spectral 
classification. 



2<r-members) have trigonometric parallaxes measured by 
Hipparcos/Tycho, these parallaxes are significant (i.e. n > 3cr ff ) 
for only 640 of these stars. Reliable Hipparcos-based distances 
to a few nearby clusters have already been determined by 
Robichon et al. ( 1999 1. But for the vast majority of the clus- 
ters under study, trigonometric parallaxes are not sufficiently 
accurate or not available at all. Therefore, we are forced to use 
indirect estimates of distances. 

The photometric approach of simultaneous determination 
of distance and interstellar extinction requires at least 3-colour 
photometry which is lacking in the ASCC-2.5. Therefore, for 
255 clusters we took the data from a list derived and newly re- 
vised by Loktin et al. (1200 II 12004 1 (LGM hereafter) who esti- 
mated distances and interstellar extinctions from homogeneous 
photometric parameters. Further, we include data of Robichon 
et al. l[T$33J - 8 clusters, Lynga fTW7\ - 19 clusters, DLAM - 
3 1 clusters, de Zeeuw et al. ( 1999 1 - 1 association, Ruprecht at 
al. dl98U - 6 associations. 

For the remaining 200 clusters we determined or re- 
vised 3 cluster distances and reddening by use of supplemen- 
tary data on spectral classes of the most probable members 
available from the ASCC-2.5 and the Tycho-2 Spectral Type 
Catalog (Wright et al. 2003 1. Further, we adopted Schmidt- 
Kaler's ( 1982 1 Zero Age Main Sequence (ZAMS) for the Main 
Sequence (MS) fitting, and A v = 3.1 x E(B - V) for com- 
puting the total interstellar extinction and colour-excess. The 
spectral class-colour/absolute magnitude calibration was based 
on Straizys J1992i . 

The references for adopted extinction and distances are pro- 
vided in the COCD for each cluster. 

In order to check the reliability of the determined distances 
and extinction, we applied our procedure to a few clusters se- 
lected arbitrary from LGM and distributed over a wide range 
of distances. The results are compared with the original LGM 
data and given in Fig. 13 We may conclude that distances and 
colour excesses derived by the spectral method coincide well 
with the published LGM data based on 3-colour photometry. 

3 For several clusters the published distance estimates do not fit the 
observed CMDs and were, therefore, recalculated (e.g. for the Pleiades 
we obtained a distance of 130 pc instead of 118 pc in Robichon et 
al. fl999l and 150 pc in LGM). 



j.z. Angular sizes 

Spatial parameters of the clusters were determined from pro- 
files of stellar density derived from star counts. For each clus- 
ter, the counts were carried out in concentric circles around the 
cluster centre. The centres were determined as the points of 
maximum surface density of the cluster members (see Sec[21>. 

As a rule, the differences between the determined coordi- 
nates of cluster centres and those listed in DLAM are small. 
Only for 21 clusters they exceed Q°.3 either due to errors in the 
cluster coordinates in old catalogs (particularly, in the Co Hinder 
list), or due to problematic definition of centres of large open 
clusters like Mel 20 (a Per), Mel 111 (Coma Berenicis), poor 
stellar groups (some objects from Platais et al. 1998 1, and 
cluster-like associations as, for example, Sco OB4. 

A general model describing the structural parameters of 
clusters was developed by King ( 1962 > and successfully ap- 
plied to globular clusters. Since a typical open cluster has a rel- 
atively small number of members and the spatial distribution 
of members is not regular, it is difficult to find a formal model 
which would, on one hand, describe structural details of a given 
cluster and, on the other hand, would be valid in general. As a 
compromise, we assumed a centrally symmetric distribution of 
the cluster members and considered only two structural com- 
ponents, a core with a radius r\ and a corona with a radius r^- 
The core radius corresponds to a distance where the decrease 
of stellar surface density stops abruptly. The corona radius (i.e. 
the actual radius of a cluster) is defined as the distance from the 
cluster centre where the surface density of stars becomes equal 
to the average density of the surrounding field (see Fig.^b)). 
For each cluster, core and corona radii were checked by visual 
inspection, and the constrains set by the corresponding VPD 
and CMD were always taken into account. Even though this 
approach is somehow subjective and time consuming, it can be 
applied to poorer clusters and thereby considerably expand the 
sample studied. 

For each cluster, the surface density was computed in con- 
centric strips of 0?05 up to 10 degrees from the cluster centre 
and for three different stellar samples: (i) all stars, (ii) cluster 
members with probabilities P > 1%, and (iii) cluster members 
with probabilities P > 61%. Nevertheless, the distribution of 
lcr members (P > 61%) was the decisive factor for the de- 
termination of the cluster radius. On average, a core radius is 
about 2.5 times smaller than a corona radius. Fig.^(b) illus- 
trates the degree of variation of F(r) with radial distance in the 
case of open cluster Stock 2. 

For 515 clusters in common, Fig. [5] compares the angular 
radii of cluster coronas determined in this paper and the cor- 
responding data compiled by DLAM 4 . Independent of the dis- 
tances, the published radii are in average lower (by about a fac- 
tor of 1.5 for "large" clusters and by a factor of 2.5 for "small" 
clusters). There can be several reasons for this bias. Often, the 
cluster sizes are empirically derived from star counts or even 
from visual inspection of cluster areas, without previous mem- 
bership determination. The area of concentration of the bright- 
est stars usually defines the cluster size, and fainter cluster 

4 The DLAM list includes a number of clusters with cluster sizes 
taken directly from the Lynga catalogue 
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Fig. 3. Comparison of derived angular radii with data compiled 
by Dias et al.( 2004 1. The solid line is the locus of equal radii. 
Crosses mark clusters with distances d > 1 kpc, whereas the 
dots are nearer clusters. The dashed line separates "large" and 
"small" clusters. 

members are lost in the rich fore- or background. In this case, 
one should speak rather on the "core radius" than on "cluster 
radius". Sometimes, a cluster size is estimated from dedicated 
observations in small sky fields. The investigations are limited 
by the detector used in the study (e.g., photographic plate or 
CCD frame), and therefore, give rather underestimated cluster 
radii. Finally, adopting the cluster sizes from the Lund cata- 
logue, one should not forget that Lynga( 1987) himself aimed 
to publish the sizes of the cluster cores. 

Our determinations of cluster sizes are not restricted by 
limited sky areas, and the information on membership is used. 
Nevertheless, the sizes derived are based on counts of relatively 
bright stars (V < 12), and thus they can be influenced by the 
mass segregation effect (see e.g. Raboud & Mermilliod 1998 
or Kharchenko et al. 2003 1. Therefore our data should be con- 
sidered as lower limits of the actual cluster sizes. 

4. Cluster kinematic parameters 

4.1. Proper motions 

For each cluster, the mean components Jl x v of common proper 
motion were computed from proper motions of the most proba- 
ble members (P > 61%). The number of such stars varies from 
3 to 178, with an average of 18 stars per cluster. For 91% of 
the clusters, the proper-motion components Ji x y are determined 
with standard errors of less than +1 mas/yr (for 42% - less than 
+0.5 mas/yr). The proper motions are derived directly in the 
Hipparcos system, for 301 clusters - for the first time (see also 
Kharchenko e. a. 125031 

Fig. |4] compares the cluster proper motions derived in this 
study with the corresponding findings of Robichon et al. il999> . 
Baumgardt et al.(2000l, Dias et al. d200ll 12002b which are 
based on the data from the Hipparcos and Tycho-2 catalogues. 
Although different techniques of member selection were ap- 
plied, the mean proper motions of the common clusters agree 



-20 20 -40 -20 20 -40 -20 2( 
/J. x ( ref ). mas/yr fx x (ref), mas/yr /z x (ref), mas/yr 




-40 -20 



20 -40 -20 



20 -40 -20 



fj. y (ref). mas/yr 



/It y(ref). mas/yr 



Fig. 4. Comparison of cluster proper motions derived in this 
study from ASCC-2.5 data with results of (a) Robichon et 
al. fi99"9l (b) Baumgardt et al. @000l, and (c) Dias et 
al. (|2001 2002). Bars show rms-errors. The straight lines are 
the loci of equal proper motion components. The upper lines 
and left axes refer to clusters included in our sample and in the 
samples of comparison. The lower lines and right axes are for 
35 clusters common to all samples. 



quite well. Even clusters with small proper motions coincide 
within a few mas/yr. This is remarkable because in these cases 
cluster members are difficult to separate from field stars. 

4.2. Radial velocities 

It is well known that our knowledge on radial velocities (RV) 
of open clusters is much poorer than on proper motions. 
According to DLAM, RVs have been published for only 240 
of about 1700 known clusters. Moreover, these data are very 
inhomogeneous: sometimes the RV of a cluster is taken from 
measurements of only one star, sometimes the rms errors reach 
30 km/s, and some authors did not give any information on the 
accuracy at all. Fortunately, the situation will be considerably 
improved when the RAVE programme (Steinmetz 2003 ) will 
be completed in the next years, including some dedicated ob- 
servations in the Galactic plane. 

Radial velocities for only 196 clusters of our sample 
are listed in DLAM, the Lund Catalogue and Ruprecht et 
al. dl981> . In order to update the RVs of these clusters, we 
cross-identified (Kharchenko et al. 2004a I the ASCC-2.5 with 
the General Catalogue of Radial Velocities (Barbier-Brossat & 
Figon 2000). On the basis of our membership determination, 
we were able to revise the RVs for 159 clusters. Additionally, 
for 94 clusters, RVs have been determined for the first time. So, 
we can now publish RVs of 290 open clusters of our sample. 

In Fig. |5] the revised RVs for 160 clusters are compared 
with the published data collected in DLAM. The mean dif- 
ference between the "old" and "new" RVs of these clusters is 
RV r ef - RV = 0.36 + 0.88 km/s. 
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Fig. 5. Comparison of cluster radial velocities derived in this 
study with data published in Lynga (1987i, Ruprecht et al. 
( 1981 1, and DLAM. If the mean RV of a cluster is determined 
from more than two members, error bars are shown. 



5. Cluster ages 

In the preliminary version of the COCD (Kharchenko et 
al. 2003 1, cluster ages were taken from the literature. Since dif- 
ferent authors used different methods of age estimates, these 
ages represented an inhomogeneous set of data in contrast with 
the homogeneity of other cluster parameters in the catalogue of 
Kharchenko et al. (2003). In the current version of the COCD 
we implement our own isochrone -based procedure of age deter- 
mination, which provides a uniform scale of ages for all clus- 
ters. 

The contents of the COCD catalogue put several constraints 
on the theoretical models which can be used for the age esti- 
mates. Since a wide span of cluster ages is expected, the CMDs 
of many clusters in our sample should present the evolved por- 
tions of the upper Main Sequence. On the other hand, the Pre- 
MS branches, observed at relatively faint absolute magnitudes 
should be seen in the COCD cluster diagrams in rare cases of 
young and nearby clusters only. Taking the above into account, 
we have been concentrating on the Post-MS isochrones. Due 
to the rather bright limiting magnitude of the ASCC-2.5, our 
sample is biased towards local clusters with typical distances 
less than 1.5 kpc from the Sun. This means that for our sample, 
no considerable metallicity trends due to the radial gradient of 
[Fe/H] in the Galactic disk is expected, and we can limit our- 
selves to the consideration of solar metallicity isochrones only. 

5.1. Input data 

There are three recently published data sets of Post-MS 
isochrones suitable for our purposes (Lejeune & Schaerer 
I2001I called the Geneva grid hereafter; Girardi et al. 2002 or 
the Padova grid; and Yi et al. l200T1 Kim et al. l2002l and Yi et 



al. 2003 known as the Y grid). Although the isochrones of the 
Geneva grid provide the widest range of ages, they do not show 
agreement with present-day observed MSs of open clusters, 
whereas the two others do (see Grocholski & Sarajedini 2003 1. 
The T 2 -isochrones agree well with cluster MSs, and integrate 
both Pre- and Post-MS stages. Unfortunately, they are limited 
to small and moderate mass models (m S 5m Q ), and cannot 
be applied to the youngest clusters with massive and bright 
stars. Thus we are left with the Padova Post-MS isochrones, 
which show good agreement with cluster MSs for My ~ 7 
mag and are finely spaced over the age scale (Alogf = 0.05). 
The major inconvenience of the Padova grid is the relatively 
high lower limit of the age scale (log? = 6.6), which restricts 
a proper dating of young clusters. We implement the Padova 
overshooting isochrone grid with the following parameters: 
m = 0.15. .. 66 m Q ,Z = 0.019, Y = 0.273. 

The Pre-MS isochrones were derived by the use of the 



Grenoble Pre-MS tracks Internet-server of Siess et al. (|2000l. 
We computed a grid of the Pre-MS isochrones which covers the 
same scale of ages as the Post-MS grid. Their agreement with 
the Padova models is acceptable and suits our purposes. So, in 
this paper we use models with overshooting for a mass range 
O.1...7.Ora andZ = 0.02, Y = 0.28. 

Since the turn-on points are observed only in a few clusters 
of our sample, we do not use the Pre-MS isochrones for cluster 
age determination but show them in the cluster CMDs to illus- 
trate the degree of conformity of the nuclear and thermal age 
scales. Further, the Pre-MS grid was already used in this study 
for photometric selection of cluster members (see Paper I for 
details). 

5.2. The method 

We applied a simple logic for the evaluation of cluster ages: 
the age of a cluster is defined by the average age of individ- 
ual MS cluster members. In contrast to the standard method 
of isochrone fitting, this approach does not only yield objec- 
tive estimates of the cluster age but also of its uncertainty (the 
standard deviation). Moreover, the algorithm is coded easily 
and can be used in the data processing pipeline. Along with 
the iterations of membership determination, the age evaluation 
procedure can be run in parallel. This is an important feature 
because one is confronted with hundreds of clusters in a wide 
range of ages. 

The individual ages of stars are derived from their locations 
in the CMD with respect to the isochrone grid (see Fig. [6). Only 
the most probable kinematic members are considered. Further, 
in order to avoid additional uncertainties which could be intro- 
duced by age estimates of red giants (e.g., due to the treatment 
of convection, mass loss, insecure conversion to the observed 
passbands), we restrict the age evaluation to MS stars 5 . So, 
we choose a red border defined by the Termination Age Main 
Sequence (TAMS) which is shown in Fig. [(J as the thick curve 
most to the right. On the other hand, a relatively slow evolution 



5 In a few cases, however, we were forced to use Post-MS stars. 
Then we selected only those of them, which are located in the blue 
part of the Hertzsprung gap. 




Fig. 6. Illustration of the adopted algorithm for determination of 
cluster ages. In Panel (a), the thin curves are isochrones for three 
different ages (log t = 8.30; 8.60; 9.00). The dotted curve is the 
locus of unresolved binaries with mass ratio q — 1 for the age 
log t = 8.60. The thick dashed curve is the adopted ZAMS where 
the dots mark the empirical ZAMS of Schmidt-Kaler dl982l l. 
The thick curves to the right of the ZAMS are the blue (s = 0.15) 
and red (i.e., TAMS, s = 1.00) borders of the evolved MS. The 
horizontal lines indicate the EMS spread at My = 0; 1; 2. The 
age variations along these lines are shown in Panel (b) as func- 
tions of the evolutionary status s (the blue edge of the evolved 
MS is marked by the vertical line at s - 0.15). Photometric er- 
rors in (B - V) of 0.01 mag are converted to corresponding errors 
in s and shown as bars. Stars marked by letters B, C, D, E are hy- 
pothetical cluster members, where B is the brightest MS star. X 
marks the isochrone base point defined as the cross-section of the 
corresponding isochrone with the blue edge of the evolved MS. 
AMy is the magnitude range adopted for the cluster age evalua- 
tion. Whereas stars B and C will be included by the algorithm in 
the age determination, stars D and E will be rejected. 



close to the ZAMS leads to a high density of the isochrones in 
this region of the CMD and, consequently, to a low accuracy of 
the interpolation. This forces us to use only an evolutionary ad- 
vanced region of the MS where the isochrone crowding effect 
is much weaker. This "blue" border is shown as the curve to the 
right of the ZAMS and is defined below. We call the region of 
the cluster MS between the red and blue borders the Evolved 
Main Sequence (EMS). 

In order to define the blue border of the EMS, we introduce 
a parameter s as a measure of the evolutionary advancement 
of a given star (called hereafter the evolutionary status of the 
star). For a star with absolute magnitude My and colour index 
(B - V)o, the parameter s is defined as 

(B - V) - (B - V)zams 

s — > 

(B - V)tams —(B — V)zams 

where (B - V)zams an d (B - V)tams are the colour indices 
of stars of the same magnitude My on the ZAMS and TAMS, 
respectively. Since the ZAMS is not explicitly present in the 
Padova grid, we adopt the youngest isochrone of log t = 6.6 
to represent the ZAMS. At My > -3 this isochrone coincides 
well with the empirical ZAMS of Schmidt-Kaler ( 119821 1. 

In Fig. [(J) we show age profiles of MS stars for three dif- 
ferent absolute magnitudes. The profiles have been constructed 
by spline interpolations along the isochrones of different ages. 
As one can see, at small s (near the ZAMS) the profiles steepen 
considerably, making the age uncertainties due to typical pho- 
tometric errors unacceptably large. The analysis of the MS 
band over its full length (My = -8... + 5) has shown that for 
a relatively safe age determination one should consider stars 
with s > 0.15. For these stars, the ages are determined from 
the age profiles corresponding to their absolute magnitudes My 
and evolutionary status s. 

In order to define a range of absolute magnitudes AMy 
where the age evaluation seems to be reasonable for a given 
cluster, we estimate the age of the brightest member among 



Table 1. Distribution of clusters according to the number of 
members (N t ) used for age determination. Typically, the more 
stars could be used the more reliable is the age determination. 
The cases i to Hi are explained in the text. 
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the evolved MS cluster stars and consider the corresponding 
isochrone. Then, the bright limit of AMy is the absolute mag- 
nitude of this member, and the faint limit is defined as the cross- 
point between the corresponding isochrone and the blue edge 
of the evolved MS. If the brightest member falls just in the 
TAMS, then the full spread of the evolved MS is available for 
the age calculation. On the other hand, if the brightest member 
is located near the blue edge (i.e., AMy « 0), no other stars are 
included in the cluster age determination. The distribution of 
clusters over the number of the most probable kinematic mem- 
bers used for the computing of cluster ages is given in Tabled 
In Tabled there are three special cases marked by i, ii, or 
Hi where the applied procedure of age determination does not 
work. The first two groups include young clusters whereas the 
third one consists mainly of old clusters. For a remote young 
cluster, the ASCC-2.5 contains only the top of the cluster MS 
which is generally badly populated due to the IMF depletion 
in the domain of massive stars. For these clusters the Padova 
grid could not be used for our purpose, because for stars with 
My 5 -3 mag, the computed strip of the evolved MS becomes 
artificially narrow. On the other hand, different effects like vari- 
able extinction, stellar binarity, rotation etc. increase the ob- 
served spread of stars around the MS. Therefore, even small 



photometric effects could lead to a situation that no members or 
remote young clusters would be present in the narrowed EMS 
area. In total, 53 of the young clusters do not fulfill the criteria 
set by the procedure. As a compromise, we consider the bright- 
est proper motion members approaching the EMS-area either 
from the blue (group i in Tabled or from the red (ii) direction 
and treated them in the procedure of age determination as they 
were just located at the EMS borders (usually, a shift of less 
than 0.05 mag parallel to the (B - V)o axis is considered accept- 
able). This means that the cluster ages can be overestimated for 
the group i or underestimated for the group ii. Wherever it is 
possible, we prefer the red-edge ages as less affected by ran- 
dom errors. 

For 14 old clusters of the third group (Hi) we found that 
their turn-off points are fainter than the limiting magnitude of 
the ASCC-2.5. Since only sub-giants or red giants in these 
clusters are present in the ASCC-2.5, we are not able to deter- 
mine their age, so we accept, instead, the published values. For 
1 1 other clusters, the turn-off points are still seen at the faint 
ASCC-2.5 limit but due to large photometric errors at these 
magnitudes, the observed CDMs are very vague. Therefore, 
ages were estimated manually by isochrone fitting. Both cases 
fall into the group (Hi). The remarks on these special cases 
and decisions made in estimating cluster ages are given in the 
notes file added to the COCD. In order to illustrate the con- 
formance of the derived ages with the structures of observed 
CMDs, the corresponding isochrones are shown in the CMDs 
of the Cluster Diagram Atlas. 

In total, the ages were determined for 506 out of 520 clus- 
ters of our sample, and for 196 of them age estimates are given 
for the first time. 

5.3. Accuracy of the results and comparison with 
another scale of cluster ages 

Although the proposed method of individual age determina- 
tion is sufficiently flexible and can be applied to poorly pop- 
ulated and sparse CMDs, there are several sources of uncer- 
tainty which could considerably affect the results. We can di- 
vide them into random errors and biases. Random errors create 
a scattering in a cluster CMD and stem from random photo- 
metric errors, variable extinction, stellar spots etc. Some biases 
like stellar rotation or unresolved multiplicity have strong ran- 
dom parts like the orientation of rotation axes, or the actual 
distribution of the component masses. They produce a quasi- 
random scattering of stars in the CMDs. The biases due to the 
adopted cluster distance and average reddening affect system- 
atically the age estimates of one particular cluster but they act 
more or less randomly if we consider a large sample of clusters. 
Finally, there are purely systematic effects which influence the 
age estimates of all clusters in the same way, e.g. uncertainties 
of the model grid which arise basically from the underlying 
physics or from neglecting particular evolutionary phases. 

The uncertainty of a random or quasi-random impact on 
the cluster age determination can be estimated from the age- 
absolute magnitude relation (AMR) at the red edge of the EMS. 
Since the general shape of the AMR is similar at the red and 



blue edges, one could, in principle, consider the blue AMR 
edge, too. But due to technical reasons like an insufficient ac- 
curacy of published input data and slow evolution near the 
ZAMS, the red edge of the AMR suits better to our purposes. 
From rms errors in absolute magnitudes cr Mv , the age uncer- 
tainty crjog, can be evaluated as crj* - Y ' °% v ' ^ ere 7 15 
the AMR slope y = dlogt/dMy- For the Padova isochrones, 
this slope varies from 0. 1 to 0.4 within the complete MS range 
(My = -8, +4 mag). In further estimates we use an averaged 
slope of yMs = 0.26 + 0.10. 

Let us first estimate the impact of uncertainties in cluster 
distances onto the derived cluster ages. The effects we con- 
sider include uncertainties in distance modulus and reddening. 
Generally, the clusters of our sample are located within 1 .5 kpc 
from the Sun, i.e. at galactocentric distances where no sub- 
stantial radial metallicity gradient has been detected (see e.g. 
Andrievsky et al. 2002). Thus, we can adopt the solar abun- 
dance of heavy elements in clusters of our sample and neglect 
metallicity corrections. A spread of cluster metallicity could 
only arise from a Galactic disk inhomogeneity which according 
to Vereshchagin & Piskunov ( 1992 ), is about A[Fe/H] a 0.1, or 
o~m v ~ 0.1 mag if converted to absolute magnitudes. Further, 
a typical accuracy of about 0.2 mag is expected for a cluster 
distance modulus (Subramaniam & Sagar 1999 1, and a redden- 
ing uncertainty of a few hundredths of a magnitude is derived 
from the LGM data. Then, the total rms error due to uncer- 
tainties in metallicity, interstellar extinction and distance mod- 
ulus does not exceed 0.3 mag. The corresponding age error is 
o- togt *0.08. 

The unresolved multiplicity of MS stars is another effect 
which may influence the age determination. For close unre- 
solved binaries (UBs) an offset of up to AMy = 0.75 mag ex- 
ists compared to single stars. At first glance, the displacement 
produced by an unresolved component introduces a systematic 
error in the cluster age determination because it always shifts 
the UB up and to the right in the CMD. For the adopted tech- 
nique of averaging individual ages it could be regarded as a ran- 
dom effect, however, with no strong bias component. Indeed, in 
the vicinity of the MS turn-off/turn-on points, the loci of sin- 
gle stars and UBs of the same age do cross. This produces a 
symmetry in the UBs spread around the isochrone of a single 
star (see Fig. [(J). In the plot we compare two isochrones cor- 
responding to populations of single stars and UBs in a cluster 
with an age of logf = 8.6. To show the effect at its largest 
amplitude we use the extreme case of components of equal 
mass. The resulting symmetry of the UB locus with respect 
to the crossing-point (B) in Fig. [6^ is evident. One can see, 
that unresolved binaries being treated as single stars produce 
ages both lower than the cluster age (those residing above point 
B) and higher than the true cluster age (below point B). Since 
we consider in our approach only the range of evolved stars 
(those brighter than My(X)), the mentioned symmetry is not 
violated by fainter unevolved UBs. Thus, no systematic age 
bias is expected from the unresolved binaries effect in our dat- 
ing technique. On the contrary, for the isochrone-based dating 
this effect carries a potential danger: ages are underestimated, 
when evolved bright binaries are treated as single stars. Since 
AMy = 0.75 is the maximum magnitude spread produced by 
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Fig. 7. Comparison of the derived cluster ages with data pub- 
lished in LGM. Dots indicate clusters where only one star 
(N t — 1) is used for the age determination, whereas crosses and 
filled circles are for clusters with N, = 2 and N, > 2, respec- 
tively. The circle sizes are scaled by the number of stars used 
for the age determination of a given cluster (cf. Table 0. For 
clusters with N, > 2, the bars show rras-errors of the averaged 
ages derived in this paper. For clusters with N, = 2, the bars in- 
dicate the scatter and are shown without hats. Triangles are the 
clusters of group (f), whereas the upside-down triangles mark 
the group (if). The dashed line is the locus of equal ages. The 
solid line shows a regression solution obtained for the complete 
sample. 

UBs in the vicinity of the single-star turn-off point, we estimate 
that this effect introduces cr log/ < 0.19 to the average cluster 
age. 

We do not consider the impact of mass loss onto the age de- 
termination since this effect is implicitly taken into account by 
the isochrones involved. Although stellar rotation can change 
the observed absolute magnitude of the MS stars, the effect 
is less important than the one introduced by unresolved mul- 
tiples. According to Maeder & Peyntremann ( 1970 1, for criti- 
cally rotating stars (near the break-up velocity) with masses of 
1.4-5 ot , the offset in brightness is less than +0.6 mag (i.e. 
cn og , < 0.16). This effect decreases rapidly with slowing down 
the rotation, and at a rotation velocity of 0.8 of the critical one 
the offset is \AM V \ $ 0.2 (cr log( < 0.05). Due to the random 
orientation of the aspect angle it has a random character. 

Also, the uncertainties of the age calibration due to the 
hooks of the isochrones near the red border of the EMS are 
small, since the time of the overall contraction after the hydro- 
gen core exhaustion is relatively short. According to the sim- 
ulations which take into account this evolutionary phase, log t 
of an individual star within the hook area changes by less than 
0.02. 



We conclude that the accuracy of age determination tor 
clusters in our sample is about o"i oer = 0.20 - 0.25 and com- 
pare the results with uniform data on cluster ages recently pub- 
lished by LGM. The cluster ages in LGM are derived via the 
method of isochrone fitting and they have the same theoretical 
basis (i.e. the Padova grid). The authors make use of photoelec- 
tric UBV observations without an explicit selection of cluster 
members. Totally, the LGM and our samples have 255 clus- 
ters in common, with 52 of them in the groups (f) and (ff) (cf. 
Table [0- Since both methods use the same isochrone grid, the 
ages of a single cluster differ only due to the applied algorithms 
and the selected observations. 

The results of the comparison are presented in Fig. As 
both coordinates are subject to rms errors, a simple linear re- 
gression analysis is inappropriate. A least square bisector solu- 
tion as proposedby Isobe et al. ( 1990 1 should rather be adopted. 
For all clusters in common, this solution of the equation 

logr = c + c! • log ftcM, 

yields (c , d)=(-0.48 + 0. 18, 1 .07 + 0.02). It is shown in Fig.|7] 
as a solid line. 

The mean age difference turns out to be Alogf=-0.06 + 
0.02, which is small compared to the standard deviation of 
o"Aiogr=0.32. According to Fig. [7] a relatively large spread is 
caused by 27 outstanding clusters above the regression line and 
4 clusters below this line. Checking these clusters in detail, we 
found that only for a few very distant clusters the large de- 
viations can be explained by the low accuracy of the cluster 
CMDs near the limiting magnitude of the ASCC-2.5. For the 
majority, the different stellar content (i.e. adopted members) 
between LGM and this paper is the decisive factor. In case of 
clusters above the regression line, bright stars in the cluster ar- 
eas were considered by LGM as cluster members, but accord- 
ing to their kinematics we found that they are non-members. 
Therefore, we may assume an underestimation of the true ages 
for these clusters by LGM rather than an overestimation by us. 
In the case of the clusters below the regression line, the ef- 
fect is reversed: bright stars within the central cluster area were 
rejected by LGM as non-members, although they fulfill all se- 
lection criteria we adopt for members. 

Excluding these 31 clusters from consideration, we ob- 
tain for the mean age difference and its standard deviation 
Alogf=0.01 + 0.01 and <TAi ogf =0.20, respectively. We may 
conclude that the applied technique allows to derive age es- 
timates with an accuracy which is comparable with the accu- 
racy of the classical method of isochrone fitting. According to 
the discussion above, the small but significant systematic bias 
(ci = 1.07) can be explained by a slight underestimation of 
the LGM isochronic ages of older (log t > 8.3) clusters due to 
the UB effect. On the other hand, since the lower limit of the 
Padova isochrones is set at logf = 6.6, this yields a bias by 
somewhat "compressing" the left side of the above relation and 
might increase cj, too. 

6. Conclusions 

Starting from a homogeneous sky-survey, we determined a 
number of astrophysical parameters of Galactic open clus- 



ters. 1ms sky survey, i.e. the ASCC-2.5, allowed us to per- 
form an unbiased rediscussion of the membership in clusters 
(Kharchenko et al. I2004b> . a necessary requirement for the 
characterisation of the clusters. The results of this paper are 
published in the form of an open cluster catalogue, which con- 
tains for all of the 520 clusters investigated: heliocentric dis- 
tance, interstellar extinction along the line of sight, angular 
size, mean proper motion and radial velocity, and age. 

For 200 clusters heliocentric distances have been newly de- 
termined via MS -fitting. We calibrated our results by a com- 
parison with a subset of clusters in Loktin's (Loktin et al. 2001 
2004 1 sample, distributed over a wide range of distances. 

For 301 clusters mean proper motions in the Hipparcos sys- 
tem have been derived for the first time. They are based on 
the individual proper motions of the most probable members. 
For the remaining clusters proper motions are compared with 
the results of Robichon et al. I ll 9991 1. Baumgardt et al. (FJOOOl, 
and Dias et al. ( 20011 l2002t and good agreement within a few 
mas/year has been found. 

For 290 clusters in our sample, mean radial velocities could 
be determined based on identifications of our members in the 
catalogue of Barbier-Brossat & Figon ( 2000 1, or retrieved from 
the literature. For 94 clusters radial velocities are determined 
for the first time, for others our findings were compared with 
Lynga d!987i . Ruprecht et al. (1981 1, and DLAM, and good 
agreement was found. On the other hand, there are 230 clus- 
ters in our sample, for which not even a single radial velocity 
measurement of any member is available so far. 

Considerable effort has been put into the determination of 
cluster ages. In the end, ages could be determined for 506 out of 
520 clusters, of which 196 are first estimates. We compared our 
determinations with the work of Loktin et al.( 2001 2004 1 and 
found good agreement. For 31 cases with large discrepancies 
between this paper and LGM, these can be explained by the 
different membership criteria between the two papers. 

Angular sizes of the cluster cores and coronas have been 
newly determined. This is the area where the full-sky cover- 
age of the ASCC-2.5 is very helpful. It makes the size deter- 
mination free from selection effects, such as limited field, star 
counts or visual inspection of the sky without taking into ac- 
count membership criteria, and others. In this paper we pro- 
ceeded the following way: from our new membership list, the 
new centre and the radial density distribution function were de- 
termined. Then by visual inspection of each of the clusters, 
core and coronas radii were determined. Despite of the rel- 
atively bright limiting magnitude of the sky survey we used 
(ASCC-2.5 with a completeness limit at V — 11.5 and a limit- 
ing magnitude at V = 14.0), we find that the angular sizes of the 
cluster coronas are systematically larger than in earlier determi- 
nations in the literature. If mass segregation in open clusters is 
an important issue (see e.g. de Grijs et al. 2002 1, our determi- 
nations of the cluster coronas are probably only lower limits of 
the true sizes. 

To improve this situation, a homogeneous sky survey with 
fainter limiting magnitude is needed. Although such surveys 
are available in infrared photometry (e.g. 2MASS), a survey 
is missing in the optical and ultraviolet regimes, it is missing 
as far as proper motions are concerned, and it certainly can- 



not be expected soon tor radial velocities. More than a decade 
from now, the ESA project GAIA is supposed to provide a sur- 
vey which will fulfill all the requirements stated above. In the 
meantime, progress has to be made in the field of ground-based 
astronomy in order to derive a homogeneous, bias-free sur- 
vey in multicolour photometry (SEGUE in the Galactic plane), 
in proper motions (beyond UCAC), and in radial velocities 
(RAVE). 
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Appendix A: The Catalogue of Open Cluster Data 

The Catalogue of Open Cluster Data (COCD) exists in machine 
readable form only and can be retrieved from the CDS online 
archive. 6 The catalogue consists of three files: format descrip- 
tion (ReadMe), the main table with the derived parameters and 
literature data, and a notes file. The necessary details on the 
data format can be found in the ReadMe file. In order to in- 
form the reader on the data scope included in the catalogue, we 
describe here the main table. The table contains data on 520 
clusters and consists of 520 lines, and 30 columns. The column 
description is given in Table IA.fl 

Appendix B: The Open Cluster Diagrams Atlas 

The Atlas (OCDA) presents visual information on the data used 
for the determination of cluster parameters, on the quality of 
member selection, and on the accuracy of the derived param- 
eters. The OCDA consists of 520 PostScript plots stored as 
gzipped files (i.e. one file per cluster), which will be available 
in electronic form only via the CDS. In order to get an easy ac- 
cess to the catalogue data, the file name includes the sequential 
number of the corresponding cluster in COCD and the cluster 
name (column 2 in the main table). Each plot contains a header 
and five diagrams used in the reduction pipeline. An example 
is given in Fig. IB. li fer the Coma Berenicis cluster. 

In the header we provide the equatorial coordinates of the 
newly determined cluster center, the cluster number in the 
COCD and the most common designation of the cluster. The 
panels present spatial, kinematic and evolutionary information. 
The upper row: the left panel is a sky map of the cluster neigh- 
bourhood constructed with stars from the ASCC-2.5, and the 
right panel is the colour-magnitude diagram. The bottom row: 
the left panel shows radial profiles of the projected density, the 
middle panel is the vector point diagram of proper motions, and 
the two right panels are "magnitude equation" (ji xy -V relation) 
diagrams. 

The sky map: A blue cross marks the adopted cluster centre. 
Pluses are centres of all clusters located within the displayed 

6 ftp : cdsarc . u-strasbg . f r/pub/cats , 
http : //vizier .u-strasbg . fr 



Fig. B.l. An example of the atlas: A plot for the open cluster Coma Berenicis. 



area and taken from the literature, i.e. the large plus in blue is 
for the given cluster, the smaller magenta ones (absent in the 
example) are for other clusters. The large circles are the bor- 
ders of the cluster core (solid curve) and of the corona (dashed 
curve). The small circles are stars; their size indicates (only in 
this panel) stellar magnitude. The bold circles are Icr cluster 
members (P c > 61%), i.e. black circles are the members in 
the core area, red circles - members in the corona. The stars 
located outside the corona are displayed in cyan. 

The colour-magnitude diagram: In this and the following 
diagrams stars are marked as coloured dots, lcr-members are 
shown with error bars. Stars used for the age determination are 
marked as bold magenta circles. Colours of the curves: ma- 
genta for the adopted Post-MS and Pre-MS isochrones, cyan 
for the borders of the Evolved Main Sequence and the "mean" 
location of the Hertzsprung gap, blue for the borders used in 
the photometric selection (see Paper I for details). Additionally, 
cluster parameters selected from the COCD are plotted where 
the rms errors are given in parentheses, and the references to 
the published data in the brackets. The symbol (#) in the line 
with the average age (log t) indicates the number of stars used 
for the age determination. 

The density profiles diagram: The curves show the distribu- 
tions of stars with angular distance from the cluster center. The 
green curve shows all stars in the cluster area, the magenta and 



black curves are for 3<x- and Icr- members, respectively. The 
adopted core and corona radii are shown by solid and dashed 
lines. 

The vector point diagram shows the proper motion distri- 
bution of ASCC-2.5 stars in this sky area. The error bars are 
shown for the cluster members. 

The magnitude equation diagrams show the distribution of 
the proper motions versus magnitudes. The horizontal lines 
correspond to the average proper motion of the cluster, the bars 
indicate the rms errors of proper motions for the cluster mem- 
bers. 
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lable A.l. Contents or the COLD main table 



Col Label Units Explanations 



1 No 




Sequential number 


2 Name 




NGC, IC or other common 






designation 


3 RA 


h 


RA J2000.0 of the center 


4 Dec 


deg 


Dec J2000.0 of the center 


5 1 


deg 


Galactic longitude ol the center 


6 b 


deg 


Galactic latitude of the center 


7 Rco 


deg 


Angular radius of the core 


8 Rcl 


deg 


Angular radius of the cluster 


9 RV 


km/s 


Average radial velocity 


11) eKV 


km/s 


Error in RV 


11 nRV 




Number of stars used for RV 






calculation 


Iz rJVIX 


mas/yr 


[i a cos S: average proper 






motion in RA 


13 erMx 


mas/yr 


Error in PMx 


14 rMy 


mas/yr 


lis' average proper 






motion in Dec 


15 eFMy 


mas/yr 


Error in PMy 


16 PM1 


mas/yr 


jU ( cos b: average proper 






motion in 1 


17 PMb 


mas/yr 


Hi,: average proper motion in b 


18 Nls 




Number of most probable (lcr) 






members 


19 d 


pc 


Distance from the Sun 


20 E(B-V) 


mag 


Reddening 


21 V-Mv 


mag 


Apparent distance modulus 


22 source(d) 


— 


Source of distance and E(B - V) 


23 logt 


logyr 


Logarithm of average age 


24 Nt 




Number of stars used 






for the calculation of logt 


25 RVref 


km - 


RV from literature 


26 eRVref 


km/s 


Error of RVref 


27 source(RVref) 




Source of RVref 


28 logtref 


logyr 


log t from literature 


29 source(logtref) 




Source of logtref 
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